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I. 


INTRODUCTION 


Previous  work  on  the  problem  of  diffraction  bv  a  nlane  annular 
sector  (see  Fiqure  11  has  heen  with  the  exact  eiqenfunction  solution 
[1,?]  hut  results  have  heen  presented  for  only  the  quarter  plane  ii.e., 
k2  =  1/2).  The  eigenfunction  solution,  of  course,  requires  a 
larqe  number  of  terms  to  ohtain  a  converqed  result  and  this,  as  well 
as  the  complexity  of  the  solution,  restricts  its  usefullness.  As 
stated  bv  Satterwhite  [?],  "the  solution  is  qeneral  and  needs  some 
simplification  before  it  can  become  a  useful  tool  for  diffraction 
studies. " 


The  ourpose  of  this  paner  is  to  make  a  simplification  in  the 
exact  eigenfunction  solution  so  that  the  solution  may  be  more  easily 
used.  Not  only  are  the  eigenvalues  found  for  the  quarter  plane  as 
has  been  done  previously,  but  the  first  few  contributing  eioenvalues 


1 


are  found  for  almost  any  anqular  sector.  These  results  are  presented 
in  oraohical  form  for  easy  use. 


The  second  purpose  of  this  paper  is  to  combine  the  eigenfunction 
solution  with  UTD  wedqe  diffraction  theory  to  obtain  the  current  on 
the  anqular  sector  usinq  only  two  or  three  eiqenfunctions.  The  use 
of  only  a  few  eiqenfunctions  is  possible  because  the  eiqenfunction 
solution  is  used  to  find  the  current  only  in  the  reqion  very  close 
H.e.,  r=C.lx)  to  the  vertex.  The  UTD  is  then  used  to  determine  the 
current  everywhere  else.  By  contrast,  Satterwhite  [?]  used  fifty 
eiqenfunctions  to  calculate  the  current  out  to  only  one  wavelength 
from  the  vertex. 

Results  will  be  shown  for  current  distributions  on  angular  sectors 
of  various  shapes  (i.e.,  O.Uk  <0.91  and  results  will  also  be  shown 
for  the  far  field  when  the  anqular  sector  is  illuminated  by  a  short 
dipole.  The  angular  sector  will  then  he  used  alonq  with  superposition 
and  the  UTD  to  find  the  current  on  four-sided  thin  plates. 

II.  the  SPHERO-CONAL  COORDINATE  SYSTEM 

The  sohpro-conal  coordinate  system  is  one  of  the  few  coordinate 
svstems  in  which  the  vector  wave  equation  is  separable  and  it  is  in 
this  svstem  that  the  plane  annular  sector  is  a  coordinate  surface. 

The  sphero-conal  system  is  described  by  three  coordinate  surfaces 
which  are  a  sphere,  an  elliptic  cone  and  an  elliptic  half-cone  as 
illustrated  in  Figure  2.  A  plane  angular  sector  is  the  e=-rr  coordinate 
surface  of  the  sphero-conal  system  and  its  angle  is  determined  by  the 
ellipticity  parameter  k^  of  the  system.  The  r,e,<f>  coordinates  are 
related  with  the  cartesian  x,y,z  by  the  following  equations: 


x  =  rcosejl  -  k’^cos^ 


v  =  rsinf)sin<}> 


m 

lop 

z  -  rcoS(j)l/1-k  cos  0 

w Here  k*'7  =  l-k?,  0  £  k?  £  1,  0  <_  9  <n  and  0  <_  <p  <_  ?ir.  The  9  and  i> 
anoles  are  not  the  same  as  the  anqles  of  the  spherica1  coordinate 
system.  A  ooint  on  the  angular  sector  (Figure  l)  has  coordinates 
rs,es,<t>s  where  9$  =  n  and  (t>s  =  cos”^  (-^4--).  The  element  ds  of  the 
surface  on  the  angular  sector  is  given  by: 

k1  sin<l> 

ds  =  r  dr  - d(J>  .  (’1 

J  1-k  ,'>cos'><t»s 

The  surface  0=it  and  d>=0  as  well  as  the  surface  0=tt  and  <b=ir  result 
in  half  planes.  For  a  more  detailed  discussion  of  the  sphero-conal 
system  see  Satterwhite  [?  ]  and  Kraus  and  Levine  [3,4]. 

III.  EXAfT  SOLUTION  FOR  FIELOS  AND  THE 
ANGULAR  SECTOR  CURRENT  DISTRIBUTION 

The  IT  field  due  to  the  current  density  71  is  determined  by  the 
dvadic.  Green's  function  which  we  must  determine.  The  IT  field  expression 
is  of  the  form 

E(R)  =  jutf  /  r(R,R').J(R'ldv  .  (31 

v 

R”  is  a  field  point,  IT'  is  a  source  point  and  v  is  the  source  region. 

In  Satterwhite 's  paper  [1]  we  can  see  that  the  Green's  function  is 
given  by 
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where:  k-  is  the  wave  number 

g  is  the  symbolic  index  which  means  the  class  of 
the  odd  or  even  solution  of  the  problem 

A„  is  the  normalization  constant, 
q 


The  vector  wave  functions  are 
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The  Z v (icr )  function  is  the  spherical  Bessel  Function  where  the  super¬ 
script  TI  means  that 


7y(.er>  = 


and  the  suDerscript  I  means  that 

V<r)  =  Jv(icr)  .  (9) 

The  0(a)  and  $(*)  functions  are  Lame  functions  defined  by  the  scalar 
wave  equations 

J  l-k2C0S2A  ^  (J  l-k2cos?A  +[(v(v+l)k?sin?A+y]e  =  0  (ini 

and 

Tl-k’^cos'><t>  ( JT-k ’  ?cos?<i>  +  [(  +1  )k l?sin?'t)-u It*  =  n.  (ill 

The  two  Equations  MO)  and  Ml)  are  coupled  throuqh  the  eiaenvalues 
v  and  y.  The  H  field  is  qiven  by  the  curl  of  the  IT  field  as  follows: 

AM  =  —  VxFffH  =  -V  X  r  TOy’l.TNR'ldv  (121 

ojy  \ 

and  the  current  of  the  anqular  sector  is 

^J(IT)  =  nx  tJRf  r,  fl=TT,0«t><n)  -  H(r ,  fi=  it, ^  .  (131 

IV.  EIGENVALUES  AND  THE  0(A)  AND  4>(<t>)  FUNCTIONS 

From  the  above  expressions  we  can  see  that  we  always  need  the 
first  even  and  odd  0(A)  and  $(it>)  functions  and  the  eigenvalues  associ¬ 
ated  with  them.  Ry  using  the  infinite  continued  fraction  method  as 
in  [l,?],  the  first  three  pairs  of  the  main  contributinq  eigenvalues 
u  and  v  were  found.  These  are  plotted  in  Fiqure  3  where  the  ordinate 
is  k  and  tbp  abscissa  is  either  (v  or  ji)  for  0<v<_1.5. 

The  qeneral  expressions  for  the  G(&)  and  functions  are 

O^fA)  =  Z  Am  cos(?m  -  ^)A  (m=. ,0, 1 1  (14) 
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(m-. . ,0, 1 . . . ) 


0Q(  ft)  =  y.  Am  sin(?m  -  ^)  ft 
<j>ef  «+>)  =  y  Bm  cos  m<t>  fm=0,l  ,?,3, . . . ) 

*0U)  =  T-  Rm  sin  m*  (m=l,?,1,...) 


ns  i 


(16) 


(17' 


Knowinq  the  eiqenvalues  p  and  v,  the  coefficients  and  m av  he 

P  mm 

found  for  0.1<_k  <0.9  by  the  continued  fraction  method  used  in  [  l,7] 

for  the  auarter  plane  (k  =1/?)  case.  The  procedure  is  tedious  and 

laborious.  The  results  are  presented  in  Fiqures  The  coefficients 

A  and  R  mav  be  determined  from  the  fiqures  with  sufficient  accuracy 
m  m 

to  use  in  Equations  and  Thus*  with  the  aid  of  the  fiaures, 

the  current  in  a  reoion  close  to  the  vertex  mav  rather  easily  he 
determined  for  almost  any  anqular  sector  as  we  shall  see  in  the  followina 
sections. 


V.  E  FAR  ZONE  FIELD  APPROXIMATION  FOR  UNIT 
DIPOLE  SOURCE  LOCATED  CLOSE  TO  THE  TIP 


An  investigation  of  the  eigenvalues  v  and  u  shows 

p 

estinq  results.  First  we  note  that  for  any  k 


some  inter- 
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These  three  eiqenvalues  v]0p  anrl  ]  have  values 

from  0  to  1.5.  If  we  examine  Equations  and  (7),  we 
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Figure  7.  The  B  coefficients  of  the  <>(<{>)  function  for 
the  first  odd  solution. 
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Figure  9.  The  B  coefficients  of  the  c(<j>)  function  for 
the  second  even  solution. 


<r  2v(?v+i)rrv+n.5) 
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For  small  values  of  v  the  expressions  M8)  and  (Id)  are  relatively 
larger  than  they  are  for  hiq  values  of  v.  This  means  that  in  the 
dvadic  Green's  function  the  first  few  terms,  say  three,  have  relatively 
hi  a  values  comoared  to  the  others.  So,  when  we  want  to  qet  an  approx i- 
mate  expression  for  the  far  zone  electric  field  in  the  case  of  a  unit 
dioole  source  located  close  to  the  vertex,  we  can  use  only  the  first 
few  eiaenvalues  in  the  expression  of  the  dvadic  Green's  function  and 
obtain  accurate  results.  In  Figure  10  we  see  two  locations  of  the  unit 
dipole  source.  In  Figures  11  to  15  we  see  the  far  zone  electric  field 

O 

with  0.1<k  50.9,  corresponding  to  Fig.  10a.  In  Figures  16  to  20  we  see 
the  far  zone  electric  field  for  the  second  location  of  the  unit  source 
dipole.  The  computations  in  Figures  16  to  20  are  for  angular  sectors 
with  0.15k^S0.5. 


Figure  10. 


gure  11.  The  F.e  far  zone  field  in  t! 

close  to  the  tip  of  plane 
(source  on  x  axis  and  para 


0  0.2  0.4  0.6  0.8  1.0  1.2 


Figure  1?..  The  Eg  far  zone  field  in  the  x=o  plane  due  to  a 
unit  source  dipole  close  to  the  tip  of  plane 
angular  sectors  for  k.  =0,1  to  0,9  (source  on 
x  axis  and  parallel  to  y). 


Figure  13.  The  E<t>  far  zone  field  in  the  x=o  plane  due  to 
a  unit  source  dipole  close  to  the  tip  of 
plane  angular  sectors  for  k2=0,l  to  0,9 
(source  on  x  axis  and  parallel  to  y). 
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SOUR 


i  far  zone  field  on  the  plane  angular  sectors  due  to  a  unit  sojrce 
?  close  to  the  tip  for  k*=0,1  to  0,9  (source  on  x  ax*s  ard 
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VI.  CURRENT  DISTRIBUTION  ON  THE  ANGULAR  SECTOR 


The  expressions  indicate  that  the  current  in  the  r  direction  has  the 
form  l/sin*.  This  means  that  the  current  on  the  edoes  qoes  to  infinity. 
But  we  know  the  current  distribution  due  to  geometrical  optics  and 
the  edge  diffracted  field  does  not  give  infinite  current  on  the  edqes 
when  the  incident  electric  field  is  normal  to  the  edge.  So,  we  conclude 
that  the  1/sin*  must  arise  because  of  the  tip  diffracted  field. 
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Near  the  tip  the  current  distribution  is  larqelv  determined 

by  the  SDherical  Ressel  function  .1  fieri  when  the  source  is  witMn 

v 

a  distance  rQ  >  r.  For  <r«l  the  expression  of  the  current  will  depend 
on  the  j^(icr)  which  have  small  values  of  o.  As  mentioned  earlier, 
the  small  eiqenvalues  are  the  first  two  or  three.  These  qive  the 
dominant  current  near  the  tip.  With  the  help  of  Equation  M3)  for 
the  even  eigenfunctions  the  current  distribution  is 


3el  =  B(icr)Vel 
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sm<}> 


(22) 


and  for  the  odd  eigenvalues  it  is 
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where 
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Yq  is  the  free  space  admittance.  From  the  expressions  of  the  current 
near  the  tiD  we  can  see  in  general  that  the  current  depends  on  three 
different  functions. 


In  Figures  IF,  17  and  18  we  can  see  the  form  of  those  functions 
for  an  angular  sector  with  k=0.3.  Far  from  the  tip  flcr>l)  the  tip 
diffracted  current  must  be  given  by: 


4>1  ( <t>) 

sin'(j> 


f  + 


(26) 
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Figure  16a.  The  Or-(kr)  current  distribution  for  an  angu 


Figure  16b.  The  Jf  current  distribution  as  in  Fig.  16d.  as  seen  from  a  different  aspect. 


Ficur^  18b.  The  J  current  distribution  as  in  Fig.  18a  as  seen  from  a  different  aspect. 

^  "  Y' 


whe’-e  the  functions  ♦($'  are: 


#^40  =  Ao+A^cos<t>+Apcos<>'*'+. . .  ,r>~ ' 

<J>2(4>)  =  3  ♦BjSin^^sIn?**...  (??* 


T  e  unknown  coefficients  A^  and  8^  can  be  found  bv  usina  the  EFIE 
io  generate  two  systems  of  linear  equations  using  point-matchinc.  >e 
EFIE  is 

n  x  f nC  (r)  =  7^  n  x  /  (-wV^G^  3$V  G)dS 

where  G  is  the  free  SDace  Green's  function  and  in  this  problem 
is  given  by 


"T  _ 
us 


Jtota1 


V  *°  * 


0  _<  <r<l  (see  Eo.  (BOl'j 

r 

<r  >  1  > 


'  :n' 


or  in  words  for  .r.l,  the  current  is  the  sum  of  the  geometrical 
optics  current,  the  diffracted  current  from  the  edges  and  the  tic 
diffracted  current  respectively. 

VII.  ANGULAR  SECTOR  FAR  FIELD  RESULTS 

By  using  the  EFIE  with  just  three  sample  points,  we  found  the 
three  unknowns  in  the  tip  diffracted  current  expressions  in  Equations 
(27)  and  (28).  The  integration  of  EFIE  was  terminated  20  wave¬ 
lengths  from  the  sample  point  since  the  integrated  functions  go 
quickly  to  zero  with  distance.  The  sample  points  are  typically 
located  about  10*  from  the  tip.  In  Figures  19,  20,  and  21  we  can  see 
the  total  field  for  an  infinitesimal  dipole  located  near  the  angular 
sector  and  normal  to  it.  The  distance  of  the  dipole  from  the  tip 
is  ter  =2.0  with  0Q  =n  and  <f>0  =  * /2 . 
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II.  FINITE  THIN  PLATES 


A  finite  thin  plate  has  ho  n  edges  and  vertices.  We  can  use 
the  angular  sector  solution  in  section  VI  to  obtain  the  currents  in 
the  regions  close  to  the  vertices  and  the  UTD  to  obtain  the  currents 
elsewhere.  If  the  vertices  are  far  apart  by  more  than  a  wavelength 
or  so,  we  can  neglect  the  interaction  between  the  vertices  and  use 
the  expressions  in  Sections  VI  and  VII  plus  superposition  to  obtain 
the  total  contribution  from  the  vertices.  That  is,  we  use  the  EFIE 
as  in  the  previous  section,  except  that  the  number  of  sample  or  match 
points  will  be  n  times  the  number  of  unknowns  used  in  the  angular 
sector  case. 

Results  for  a  RX  hy  9X  souare  Dla*°  me  shown  in  Fiqurp  ??. 

Hi.)  to  the  nature  of  our  solution  s»e  a  larq°  "spike"  of  current 
at  the  four  corners  in  Fioure  9?a.  The  current  in  Fioure  ?9a  mav 
comnared  wi^h  that  oresee*eri  hy  Ko  anr1  t tra  L f J  ,  and  it  will 
he  seen  that  the  results  are  s i am f icant  1  y  different,  at  the  four  corners 
since  their  results  show  no  singularity  of  the  current  at  the  corners 
for  this  size  plate.  We  believe,  that  because  the  interaction  between 
corners  is  small  (or  negligible)  for  large  plates,  it  is  inevitable 
that  a  strong  singularity  be  observed  in  current  there.  This  con¬ 
clusion  is  the  opposite  of  that  in  [5].  It  may  be  that  the  solution  in 
[5]  is  not  converged  for  large  plates. 

Figure  23  shows  the  current  on  a  3X  by  3x  plate  that  is  not  square, 
but  in  the  shape  of  a  rhombic  with  two  corners  characterized  by 
k^=0.3  and  the  other  two  by  k^=0.7.  In  both  Figs.  22  and  23  the  edge 
condition  of  the  current  is  not  shown  for  purposes  of  clarity  nor  is  it 
used  in  the  far  field  computation  in  the  next  paragraph. 

Finally,  Figure  24  shows  the  far  field  pattern  of  a  short  monopole 
or  stub  radiator  at  the  center  of  a  8x  square  ground  plane.  The  result 
is  compared  with  an  unpublished  result  by  Marhefka  who  obtained  his  re¬ 
sult  using  the  UTD  with  vertex  diffraction.  The  two  results  are  seen  to 
be  in  very  close  agreement. 
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Figure  19.  Total  electronic  field  for  five  different 
angular  sectors. 
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Figure  20.  Total  electric  field  for  five  different 
angular  sectors. 


o 


Figure  21.  Total  electric  field  for  five  different 
angular  sectors. 
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IX.  SUMMARY 


In  this  paper  we  have  presented  in  easy  to  use  graphical  form 
the  first  few  contributing  eigenvalues  for  any  angular  sector  having 
a  wide  range  of  Included  angle.  We  have  shown  that  it  is  only  these 
first  few  eigenvalues  that  contribute  to  the  current  in  a  region  very 
close  to  the  tip  Itself.  Beyond  this  region  the  GTO  may  be  used  to 
obtain  the  current  thereby  overcoming  the  need  to  use  large  numbers 
of  eigenvalues  and  eigenfunctions. 

We  then  applied  the  combined  eigenfunction-GTD  technique  to 
find  the  currents  and  resulting  radiation  fields.  To  demonstrate 
the  validity  of  the  solution,  the  radiation  pattern  of  a  monopole 
at  the  center  of  a  square  plate  Is  compared  with  that  calculated  with 
an  unpublished  GTD  solution  which  Includes  vertex  diffraction.  The 
two  independent  results  are  in  very  close  agreement. 
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short  monopole  in  the  center. 
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APPENDIX 

The  following  tables  give  the  lower  order  eigenvalues  for  angular 
vectors  with  0.1  <  <  0.9  and  for  0.0  <  v  <  1.0  and  0.5  <  v'  <  1.5. 
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MISSION 
of 

Rome  Air  Development  Center 

RAVC  plant,  and  execute*  r etearch,  development,  te*t  and 
selected  acquisition  program*  In  *u pport  of  Command,  Control 
Communication*  and  Intelligence.  (C3 1)  activities.  Technical 
and  engineering  Support  ulithin  area*  of  technical,  competence 
it  provided  to  ESV  Pn.ogA.am  Offices  (P04)  and  other  ESV 
element s.  The  principal  technical  mi&tion  area*  are 
communication*,  electromagnetic  guidance  and  control,  *ur- 
veiilance  o f  ground  and  aerotpace  object*,  intelligence  data 
collection  and  handling,  information  *y*tem  technology, 
ionotpheric  propagation,  tolid  ttate  tcience*,  microwave 
phytic*  and  electronic  reliability,  maintainability  and 
compatibility. 
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